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We consider the M = 1 Skyrme model and obtain supersymmetric skyrmion solutions numerically. 
The model necessarily contains higher derivative terms and as a result the field equation becomes 
a fourth-order differential equation. Solving the equation directly leads to runaway solutions as 
expected in higher derivative theories. We, therefore, apply the perturbation method and show that 
skyrmion solutions exist upto the second order in the coupling constant. 
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Soliton solutions in non- or supersymmetric gauge theories at large-A'^ have been playing an important role to 
understand non-perturbative phenomena in QCD. 
' A well-known example is the Skyrme model and its soliton solutions, called skyrmions 0. Witten investigated 
mesons and glue states in the large- limit of QCD by a systematic expansion in powers of and showed that 
baryons emerge as solitons with mass of 0{N) 3- the successive papers Jsj, it was shown that the resultant 
effective theory is the Skyrme model and skyrmions are interpreted as baryons jj| . 
^ ■ More recently, it was found that domain wall solutions exist in large- SUSY gluodynamics (SQCD) 5] as well as 
<^ in non-SUSY QCD These walls behave as D-branes on which the string could end if they are BPS saturated 

^\ An interesting observation is that the wall width is 0{1/N) and correspondingly heavy states are expected to emerge 
i as solitons with mass of 0{N) 0. 

■ In this context, it may be natural to consider that solitons with mass of 0{N) in SQCD could have something 
^ID \ to do with skyrmions in the supersymmetric version of the Skyrme model. The supersymmetric Skyrme model was 
\^ • constructed in Ref. 0|. The extension to supersymmetry restricts us to work on a CP{1) target space rather than 

\ . The authors concluded that the possibility of the existence of soliton solutions is not excluded when the higher 

' derivative terms are take into account. 

\ In this paper we construct supersymmetric skyrmion solutions numerically. In the supersymmetric case, the Skyrme 
^ • field equation is not second-order but fourth-order in derivatives. In general, higher derivative theories lack a lowest- 
^ ' energy state and exhibit runaway solutions along with physical solutions, which makes numerical computation unsta- 
' ble. The prescription was proposed by Simon in Ref. According to his argument, if higher-order derivative terms 
^ \ can be considered as a small perturbation, physical solutions should be Taylor-expandable around the leading order 
solution. We shall apply this perturbation method to the supersymmetric Skyrme field equation and obtain soliton 
solutions upto the second order in the coupling constant. The dependence of the skyrmion solutions on the coupling 
' constant is also examined. 

2. THE SUPERSYMMETRIC SKYRME MODEL 

In this section, we give a brief review of the TV = 1 supersymmetric Skyrme model constructed in Ref. ^^I- 
supersymmetric theories, the target manifold M must be Kahlar It was shown in Ref. ^3 that the only nontrivial 
homotopy group in four dimensional spacetime is 7r3(CP(l)) = Z. The complex projective space CP{N) is realised 
by gauging the C/(l) subgroup of SU{N), i.e. CP{N) = SU(N)/SU(N - 1) x [/(I). This amounts to replacing 
ordinary derivatives in the Lagrangian into covariant derivatives I?^ = 9^ — iV^C/ra where (i = 1,2,3) are Pauli 
matrices. The Skyrme Lagrangian is then given by 

C = Jl^iD^U^D^U) + ^tr[C/ti?,,[/, U^D.Uf . (1) 

where /tt is a pion decay constant and e is a dimensionless constant. Let us parameterise the chiral field in terms of 
the complex scalars A = {A^,A^) e with AA = A^Ai + A\A-2, = 1. Ai is related to SU{2) matrix by 
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One can parameterise the gauge field in terms of Ai as 



V,^~-[Ad,A~id,A)A] 



The Skyrme Lagrangian ^ is then written as 



- - 



where 



Now, uj = U'^d^Udx^ is an SU(2)-valued one-form and therefore the Maurer-Cartan equation holds 



which reads 



Thus, the Lagrangian Q becomes 



duj + LJ /\uj — Q , 



'-D''AD„A 



16e 
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To supersymmetrise the model, let us extend {Ai, y^) to the chiral multiplet and vector multiplet respectively 

A,^{A„ijai,F,), V^,~>{Vf„Xa,D) 



(3) 

(4) 

(5) 
(6) 
(7) 

(8) 
(9) 



where i,a = 1,2 and Fi are complex scalars, I? is a real scalar and ipai, Xa are Majorana spinors. Then, the 
supersymmetric Lagrangian is given by 



J-SUSY — 



-Df'AD^A + '-D^i^^a^J"" - JV'^'T^ad^'^V;" +FF- lAX^^,^ + i^A^V^" + D{AA - 1) 
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One can show that these are invariant under the following supersymmetric transformations: 

SA, =-e"Vaz 
Si^ai = -ia'^^e" D ^,Ai + eaFi 
SF, = -le^a^^D^^f - leM'^ 

and 

With these transformations and the field equations, the constraints AA = 1 can be extended to 

AA^l , Aipa = , AF = 0, 

and resultantly one obtains 



--{Ad^A - d^AA) + -iP'^a^^^r 



K =-iF^^+a^J''Df,A 

D = D^AD^A - ^{D^^j'^a^J'^ - V^a^ad^'^^") . 



(10) 



(11) 
(12) 

(13) 



(14) 
(15) 
(16) 

(17) 

(18) 
(19) 
(20) 
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Setting ipa = Fi = 0, one arrives at the bosonic sector of Eq. H10|l . 

1 



jOsusy = -^D^^AD^A+^ 
8 8e^ 



(21) 



The last term is fourth-order in derivatives. There are other possible fourth-order terms we can add to the Lagrangian, 
which turn out llOl 



UAUA - {D^'AD^Af 
The most general supersymmetric Skyrme model is thus given by 



^SUSY — 



D^'ADu,A 



1 

8^ 



if^^ + {D'^AD^Af'^ + /? { UAUA - {D^'AD^Af^ 



(22) 



(23) 



Although we have simply set Fi = to get Eq. (|23|l . the Lagrangian (|10|l contains derivatives of Fi and hence Fi 
becomes a dynamical field. However, since the dynamical solution does not cancel the higher derivative terms to 
recover Eq. ijSJ), the Lagrangian would be still given in the form of (|23|) . 

We consider spherically symmetric solutions with the topological charge 1. Let us impose the hedgehog ansatz 



U = cos /(r) + IT ■ n sin /(r) 
where n = x/r. In terms of A^ it corresponds to 

Ai = cos /(r) + i cos 6 sin /(r) , A2 — ie^"^ sin 9 sin /(r) . 
The topological charge in the Skyrme model is defined by 



247r2 



(24) 
(25) 

(26) 



where we have used the relation Q in the equality. Inserting the ansatz H25(l into H26(l and using the boundary 
conditions /(O) = tt, f{oo) = 0, one obtains Q = \. 

We insert the ansatz (|25|) into the Lagrangian H23|) . Then the static energy is given by 



dx ; 



1 



l\2 



2 sin^ / 



^ I 2f sin 2/ 



(27) 



12 ' x2 J 15 ' a;2 J 12 

where we have introduced the dimcnsionless variable x = ef^rr and the prime denotes a derivative with respect to x 



3. FIELD EQUATIONS 

The field equation can be obtained by taking a variation with respect to f{x) 

4(a + (3) 



x^f" -2xf' + sm2f 



+f3 



a;2/(4) + 4a;/(3) _ 4/" cos 2/ + 4(/')^ sin 2/ 



2/" sin^ / - Gx^fff" Mf'f + iff sin 2/ + 
4 sin^ / sin 2/ 



sin / sin 2/ 



0. 



(28) 



This equation contains third- and fourth-order derivative terms. In the theories which contain higher-derivative terms, 
there exists no lowest-energy state no matter how small their coefficients are. Thus, if one solves them directly, one 
would end up wi th p icking up unphysical runaway solutions. The perturbative method to avoid this problem was 
proposed in Ref. jllj and has been applied for theories such as gravity and knotted solitons [l3| . 

Let us apply the perturbation to solve the equation (|28|) . The Skyrme model is a truncated theory of derivative 
expansion. Therefore, the higher-derivative terms would be treated as small corrections. Let us expand the Skyrme 
field in the order of the coefficient of the higher-derivative terms /?, 



L 

E 



/(x) = /o(x) + >>'/z(x)-fO(/3^+i) 



(29) 
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We compute solutions upto the second-order in i.e. L = 2. 
From Eq. (|28() . the zeroth-order equation is given by 



0(/3°) : h^f'^ - 2r/^ + sin2/o + ^ \-Ar{J'^f + {f'^f sin2/o + 

O T 



sir? fo sin 2/o 



where 



The first-order equation is given by 



where 



S3 



Oi/3) : /^o/r - + S2/1 + S3 + = 

5 



= 2r + ^ ( 6rV^/^' + Or/^^ - sin 2/o) 

2/0 + ^ ( 2/^' sin 2/0 + 2/^2 cos 2/0 + wo) 

4 sin^ fo sin 2/o 



S2 = 2 cos 



S4 = 2(sin^ fo 



+ 4r/('^ - 4/^' cos 2/0 4/^2 sin 2/o 

3f /o )/o - 4r/o + /o sm2/o -f ^ 



and 



Wo 



= ^{sin^ 2/0 + 2 sin^ /„ cos 2/o) . 



The second-order equation is given by 

OiP"^) : hof2 +^{2mi+ 1712) +1713 + ^7714 + 7715 =0 
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where 



mi = -6r2/o/o/2 + /osin2/o 
m2 = 2(/^sin2/o-6r/,'2^^' ' 



= 



....,u/2 + (sin 2/0/1 - Qr^f'of'^f'{ + /f/o cos2/o ~ Sr^/f /o 
'o')/2 + (2/0' cos 2/0 + i/;o)/2 + (sin 2/0 - 12r/^)/f + 4/^/{/i cos2/o 
1 . "1 



2/^2 sin 2/0 + ^ (3 sin 2/o cos 2/o - 2 sin^ fo sin 2/o) /f 



777,3 



777, 



- rVi*'^ + 4r/f ' - 4/{' cos 2/o + Sf'of'i sin 2/o + 4(2/^' sin 2/o + 2/^^ 2/o - Wo)fi 

- o^c;.2 f„ _ ,^2 f/2^ p, ^ r,, p o _ 6r f^^ _ Qr^f'J'^)f[ + (2/^' sin 2/o + 2/^^ 2/o + wo)fi 



^'Jl ^^1 ^JU ^ "J OJl 

7774 = 2(sin2/o-3r2/o')/r + 2(/^sin2/o-6r/^- 
2r/^ + 2/2 cos 2/0 -2/2 sin 2/0. 



by 



Solutions upto the second-order are then given 

/(r) = fo{x) + /3/i(x) + P^f2{x) + 0{P' 

4. NUMERICAL SOLUTIONS 



(30) 

(31) 

(32) 

(33) 
(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 
(42) 
(43) 

(44) 



We have solved the second-order differential equations 
conditions 



by the shooting method subject to the boundary 



/o(0)=7r, /o(^) = 0, /i(0) = 0, /i(^)-0, /2(0) = 0, /2(^)=0. 



(45) 



Fig- n shows the dependence of the profile function f{x) on a. The skyrmion slightly expands in size for increasing 
a. Fig. 121 shows its dependence on /?. The difference is very small but non-zero /3 contributes to expand the size of 



5 



skyrmions. Fig. O shows the dependence on a of the first- and second-order perturbed functions. The values which 
give a maximum correction to the leading order are /i ~ 1 and /2 ~ —5 for a = 0.8. If we take (3 = 0.01, each will 
give a correction O(10~^) and O(10~^) respectively. The leading order takes the value of 0(1) around the maximum 
correction values. Thus, it seems that /3 = 0.01 would be small enough to make our perturbation method valid at 
least for a > 0.8. Fig. 0] shows the a dependence of the energy density. The behavior is the same as the profile, that 
is, the skyrmion increses in size as a increases. 

The corrections gets larger as a becomes smaller, and correspondingly we have to take smaller values for (3. But 
the regime where /i and /2 take large values compared to the leading order, the derivative expansion would be broken 
down and the result should not be trusted. 

It is noted that the first-order function contributes to increase the size of the skyrmion. On the other hand, the 
second-order function contributes to decrease the size. That makes the total corrections to the leading-order solution 
even smaller. 

We have examined the dependence of the total energy on a and /3 which is shown in Fig. For increasing a, the 
energy increases monotonously. For increasing /3, the energy increases and the figure merely shifts upwords. 

In Ref. 10, it was shown that skyrmions are unstable when (3 = 0. However, for /3 7^ 0, there is a possibility that 
they are stable. Unfortunately we are not able to show if our solutions are stable since the equation for the stability 
analysis are again higher-order in derivatives. 



5. SUMMARY 



In this paper we studied the J\f = I supersymmetric Skyrme model and constructed skyrmion solutions numerically 
by the perturbation method upto the second order in the coefhcient of the higher derivative terms. The skyrmion 
depends on the coupling constants a and /3 and increases in size when these values increase. We found that the first- 
order correction contributes to increasing the size of the skyrmion while the second-order contributes to decreasing 
the size. As a result, the total correction to the leading-order solution is very small. For f3 = 0.01, the first- and 
second-order gives a correction of 1% and 0.01% respectively to the leading order solution, which should justify our 
perturbative treatment of the model. The energy of the skyrmion increases monotonously for increasing a and (3. 

It should be also noted that for /3 = 0, the supersymmetric skyrmions are unstable, but the higher-derivative terms 
could change the stability dramatically no matter how small the correction is although we have to wait for the results 
from proper stability analysis,. 

We believe that the supersymmetric Skyrme model and its soliton solutions deserves further investigation since it 
may shed light on the non-perturbative effects in large- iV supersymmetric QCD theory. In particular, quntisation of 
supersymmetric skyrmions will be necessary to be performed. 
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FIG. 1: The profile function / as a function of x for with a = 0.8, 1.0, 1.2 and /3 = 0.01 fixed. 




FIG. 2: The profile function / as a function of r for /? = 0.0, 0.1 with a = 1.0 fixed. 




FIG. 3: The perturbed profile function /i and /2 as a function of a; for a = 0.8, 1.0, 1.2. 




FIG. 5: The a dependence of the total energy with (i — 0.01 fixed. 
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